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Abstract 

On the basis of assumptions about the behavior of driver-vehicle units con- 
cerning acceleration, deceleration, overtaking, and lane-changing maneuvers, a 
gas- kinetic traffic model for uni-directional multi-lane freeways is constructed. 
Queuing effects are explicitly taken into account in an overall manner. The 
resulting model is a generalization of Paveri-Fontana's Boltzmann-like traffic 
model and allows the derivation of macroscopic traffic equations for inter- 
acting lanes, including velocity equations. The related effective macroscopic 
traffic model for the total freeway cross-section is also derived. It provides 
corrections with respect to previous traffic models, but agrees with them in 
special cases. 
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I. INTRODUCTION 



Apart from microscopic traffic models, in the last decades a number of interrelated macro- 
scopic traffic models have been proposed PHTTj. The motivations for developing these were 



to describe and understand the instabilities of traffic flow |5|-|8|,|IT|-|ilfl 



to optimize traffic flow by means of on-line speed-control systems p§-|n> 



to make short-term forecasts of traffic volumes for re-routing measures [T7|-|TT)|] , 

to calculate the average travel times, fuel consumption, and vehicle emissions in de- 
pendence of traffic volume ||J2(|, 



to predict the effects of additional roads or lanes [f2C| -|23| 



Most of these models are restricted to uni-directional freeway traffic and treat the different 
lanes of a road in an overall manner, i.e. like one lane with higher capacity and possibilities 
for overtaking. However, this kind of simplification is clearly not applicable if there is a 
disequilibrium between neighboring lanes. Therefore, some researchers carried out empirical 
investigations of the observed density oscillations between neighboring lanes or proposed 
models for their mutual influences p4|-|29|j. 

However, these are phenomenological models which treat inter-lane interactions in a 
rather heuristic way. Moreover, most of them base on the simple traffic flow model of Lighthill 
and Whitham which assumes average velocity on each lane to be in equilibrium with density. 
This assumption is not very well justified, especially for unstable traffic which is characterized 
by evolving "phantom traffic jams" or stop-and-go waves [|5|-[|,[TT]-[T3| . It is also questionable 
for lane mergings or on-ramp traffic where frequently a disequilibrium occurs. However, 
instabilities or disequilibria may decrease the freeway capacity considerably. 

An alternative approach including a phenomenological velocity equation has been pro- 
posed by Michalopoulos et al. |2£|. It bases on Payne's model []3],|]] which has been severely 
criticized for several reasons P JTI| , |30| -p5| . Therefore, we will derive a consistent macroscopic 



multi-lane model from a gas-kinetic level of description. This is related to Paveri-Fontana's 



Boltzmann-like approach (cf. Sec. II), but explicitly takes into account overtaking and lane- 
changing maneuvers. An extension to a gas-kinetic traffic model for several driving styles 
(e.g. aggressive, timid) or vehicle types (e.g. cars, trucks) is easily possible (cf. Sec. [11 C| ). 

On the way of constructing gas-kinetic traffic equations we face the problem of vehicle 
queues (platoons) forming at medium and high traffic volumes [BB]. Several solutions have 



been suggested to cope with this. Most of them base on gas-kinetic ideas, but non of them 
is fully satisfactory, yet. The approach of Andrews [p7| - |40fl mainly provides results for sta- 



tionary traffic situations. Beylich's model pl| , |42[| is so complicated that the corresponding 



macroscopic traffic equations are not any more derivable or at least not suited for efficient 
traffic simulations. For this reason, this paper does not distinguish queues of different lengths 
(like Beylich does) but only freely moving and impeded (queued) vehicles. This approach 
has some relation to the one of Lampis ||43|| , but it removes its inconsistencies. 



The resulting Boltzmann-like model allows a systematic derivation of macroscopic traffic 
equations not only for the vehicle densities on the different lanes, but also for the associated 
average velocities (cf. Sec. III). In Sec. IV it is demonstrated how effective macroscopic 
traffic equations for the total cross-section of a uni-directional road can be obtained from 
the equations for single lanes. Due to different legal regulations, the traffic dynamics on 
American freeways is different from that on European ones (which will be called "autobahns" 
in accordance with Kiihne et al. ]14|||) (cf. Sec. II). 

A summary and outlook is presented in Sec. V. 

II. BOLTZMANN-LIKE MULTI-LANE THEORY 

The first Boltzmann-like (gas-kinetic) model was proposed by Prigogine and co-workers 
^-^6[. However, Paveri- Font ana ||7[ has pointed out that this model has some peculiar 



properties. Therefore, Paveri-Fontana proposed an improved model that overcomes most of 
the short-comings of Prigogine's approach. Nevertheless, his model still treats the lanes of 
a multi-lane road in an overall manner and does not take into account queuing effects. For 
this reason, an extended Paveri-Fontana-like model will now be constructed. 

Let us assume that the motion of an individual vehicle a can be described by several 



variables like its lane i a (t), its place r a (t), its actual velocity v a (t), and its desired velocity 
v 0a (t) in dependence of time t. The phase-space density pi(r,v,v ,t) is then determined by 
the mean number Arij(r, v, Vq, t) of vehicles on lane i that are at a place between r — Ar/2 
and r + Ar/2, driving with a velocity between v — Av/2 and i> + Av/2, and having a desired 
velocity between v — Av /2 and t> + At> /2 at time t: 
Am(r,v,v ,t) 



Pi(r,v,v ,t) = 
1 



Ar Av Avo 

r+Ar/2 v+Av/2 jjo+Avo/2 

£<W (t) fdr'6(r'-r a (t)) f dv'6(v' - v a (t)) f dv' 5(v' - v° a (f)) . 



Ar Av Avn 

r-Ar/2 v-Av/2 v -Av /2 



(1) 



Here, Ar, Av, and Avq are small intervals, dij denotes the Kronecker symbol and 5(x — y) 
Dirac's delta function. The notation "/" indicates that a time-dependence only occurs in 
exceptional cases. Lane numbers % are counted in increasing order from the right-most to the 
left-most lane, but in Great Britain and Australia the other way round. (For Great Britain 
and Australia "left" and "right" must always be interchanged.) 

The phase-space density pj(r, v, Vq, t) can be splitted into a term pn(r, v, vq, t) describing 
vehicles that can move freely and a term p i2 (r,v,Vo,t) delineating impeded vehicles that 
have to move slower than desired since they are queuing behind other vehicles: 

Pi(r,v,v ,t) = pii(r,v,vo,t) + p i2 (r, v, v , t) . (2) 

Introducing the proportion Cj(r, v, vq, t) of freely moving vehicles by 

Ci(r, v, v , t)pi(r, v, v , t) = p a (r, v, v , t) , (3) 

we have the relations 

0<Ci(r,v,v ,t) < 1 (4) 

and 

Piz(r, v, v , t) = [l- c t (r, v, v , t)]pi(r, v, v , t) . (5) 
Now, we utilize the fact that, due to the conservation of the number of vehicles, the 



phase-space density pi(r,v,v ,t) on lane i obeys the continuity equation |p0|j47|j48 



+ v?{r, v, v , t) - i>i(r, v, v , t) . (6) 

The second and third term describe temporal changes of the phase-space density pi(r, v,vo,t) 
due to changes dr/dt = v of place r and due to acceleration /P, respectively. We will assume 
that the proportion Cj(r, v, Vq, t) of freely moving vehicles accelerates to their desired velocity 
v o with a certain relaxation time T so that we have the acceleration law 

f?(r,v,v ,t) = C 7 i(v -v). (7) 

The terms on the right-hand side of equation @ reflect changes of phase-space density 
Pi(r,v,vo,t) due to discontinuous changes of desired velocity vq, actual velocity v, or lane i. 
z/ 4 + (r, v , Vq, t) and uf(r,v,Vo,t) are the rates of vehicles entering and leaving the road at place 
r. They are only different from zero for merging lanes at entrances and exits respectively. 
The term 

(dpi\ Pi(r,v,t) 



■[Poi(v ;r,j£) - P i(v ;r,t)} , (8) 



where 

pi(r,v,t) = I dv pi{r,v,v Q ,t) (9) 



is a reduced phase-space density and T r m 1 s is about the reaction time, describes an adapta- 
tion of the actual distribution of desired velocities Poi(vo;r,t) to the reasonable distribution 
of desired velocities Poi(vo] r, /) without any related change of actual velocity v. 

For the reasonable distribution of desired velocities we will assume the functional depen- 
dence 



P 0i (v ; r,jt) = .1 e-N-^l 2 /^] (10) 

V 2n6 0t 



which corresponds to a normal distribution and is empirically well justified [fi9 , p5 ,20|. The 



mean value V 0i = V 0i (r,f) and variance 9 0i = 9 0i (r,f) of Poi(v ; r,f) depend on road conditions 
and speed limits. Since European autobahns usually do not have speed limits (at least in 
Germany), # i is larger for these than for American freeways. In addition, on European 



autobahns Voi increases with increasing lane number i since overtaking is only allowed on 
the left-hand lane. 
The term 

(dp,\ d 2 (, A,,/\ 

with a dimensionless, density-dependent function Ai describes a kind of "velocity diffusion" 
which takes into account individual fluctuations of velocity due to imperfect driving. For a 



detailled discussion of this term cf . Refs. [lljSO] • It provides contributions to the dynamical 
variance equations and equations for higher moments, but neither to the density nor to the 
velocity equations pOf . 

Before we specify the Boltzmann-like interaction term (dpi/ dt)- int and the lane-changing 
term (dpi/dt)\ c we will discuss some preliminaries. For reasons of simplicity we will only treat 
vehicle interactions within the same lane as direct pair interactions, i.e. in a Boltzmann-like 
manner 



50 1 . Lane-changing maneuvers of impeded vehicles that want to escape a queue 



(i.e. leave and overtake it) may depend on interactions of up to six vehicles (the envisaged 
vehicle, the vehicle directly in front of it, and up to two vehicles on both neighboring lanes 
which may prevent overtaking if they are too close). Therefore, we will treat lane-changing 
maneuvers in an overall manner by specifying overtaking probabilities and waiting times of 
lane-changing maneuvers (which corresponds to a mean- field approach, cf. Ref. [pOj] ). These 
probabilities and waiting times dependent on the vehicle densities and maybe also on other 
quantities. 



A. Overtaking 

Let, for example, P^ denote the probability that a slower vehicle could be immediately 
overtaken on the left-hand lane and Pf the analogous probability for the right-hand lane. 
In addition, let Pj = Pi (r, /) be the proportion of driver- vehicle units that desire to move 
on lane io (which corresponds to the proportion of vehicles that actually moves on lane io 
at small vehicle densities). Then, 

pf=pf(rj)= E P«(r,fl (12) 

io(>i) 
a 



is the proportion of vehicles on lane i which would prefer to overtake on or change to the 
left-hand lane i + 1, whereas 

pf=pf(rj) = £ P io (rJ) (13) 

io(<i) 

is the proportion of vehicles which would prefer to overtake on or change to the right-hand 
lane i — 1. For American freeways we will now calculate the probability pf [p~] with which 
vehicles on lane % are immediately overtaken on lane % + 1 [lane i — 1] by vehicles that also 
moved on lane i beforehand. Since -Pj + (1 — P^) [Pi~(l ~ ls the probability that only 
lane i + 1 [lane i — 1] allows overtaking, and P^P 7 ' is the probability that both neighboring 
lanes are free, we obtain 

p+ = Cj {i^(l - Pr) + [ P > + 1(1 - p> - p<)]i^ + Pr} 

= Cl [P+(l - Pr) + 1(1 +p > -p<)P?Pr] (14) 

and 

pr = Q {Pr(l - P+) + [p< + 1(1 - p> - p<)]P+P,r} 

= Cl [Pr(l - P+) + 1(1 +p< -p>)P+Pr] (15) 

for inner lanes. Here, we have taken into account that only the proportion q of freely moving 
vehicles is able to immediately overtake. Moreover, we have assumed that driver- vehicle units 
change towards their desired lane i^ if both neighboring lanes allow overtaking. If already 
being on the desired lane io = i (which is the case with probability 1 —pf —pf), they choose 
each neighboring lane with probability 1/2. If only one lane is free, this chance is taken 
irrespective of the desired lane %q. From fllil) and (|T5|) follows that the total probability of 
immediate overtaking is 

p++pr = Cl [l-(l_P+)(l-Pr)]. (16) 
Since, on outer lanes, only one lane is available for overtaking maneuvers, we have 

pf = 0, pr = Ci Pr (17) 

on the left-most lane and 



pr = 0, pf = c l P l + (18) 

on the right-most lane. 

Of course, other specifications of pf and p~ are also possible. In Europe, for example, 
overtaking on the right-hand lane is prohibited for free traffic flow (with less than 30 vehicles 
per kilometer and lane) so that 

pf = c i P l +, pr = (19) 

on inner lanes and the right-most lane, but 

p+ = 0, p~ = (20) 

on the left-most lane. However, for congested traffic (with an average velocity less than 80 
kilometers per hour) vehicles are also allowed to pass slower vehicles on the right-hand lane 
so that the situation is similar to American highways, then. 

We are now ready to specify the Boltzmann-like interaction term. For not too large 
vehicle densities it can be written in the form [20 



dpi 

at 



hit 



/ dv' dw dw W2(v,i\v' ,i';w,i')pi'(r,v' ,v ,t)pii(r,w,w ,t) (21a) 



w<v' 



~ J d v ' J dw J dwoW2(v',i f \v,v,w,i)pi{r,v,vo,t)pi(r,w,wo,t). (21b) 



Term ( |21aj) describes an increase of phase-space density pi(r, v, vq, t) by interactions of a ve- 



hicle with actual velocity v' and desired velocity v on line il with a slower vehicle with actual 
velocity w < v ' and desired velocity w causing the former vehicle to change its velocity to 
v ^ v' or its lane to i ^ i' . The frequency of such interactions is proportional to the phase- 
space density pi/(r,w,w ,t) of hindering vehicles and the phase-space density pi/(r,v',v ,t) 
of vehicles which can be affected by slower vehicles. Analogously, term ( |21b|) describes a 
decrease of phase-space density pj(r, v, v , t) by interactions of a vehicle with actual velocity 
v and desired velocity vq on line i with a slower vehicle with actual velocity w < v and 
desired velocity Wq causing the former vehicle to change its velocity to v' ^ v or its lane 
to il i. Since the interaction is assumed not to influence the desired velocities vq, Wq, the 
interaction rate W2 is independent of these. However, the interaction rate W2(v',i'\v,i',w,i) 



is proportional to the relative velocity \v — w\ of approaching vehicles. Therefore, we have 
the following relation: 

W 2 (v', i'\v, i; w, i) = pf\v - w\5 ilii+1 )5(v' - v) (22a) 
+ Pi\v-w\5 i >(i^ 1 )5(v'-v) (22b) 
+ {1 - pi)\v - w\5 n 5{v' - w) . (22c) 

Term (|22a|) describes an undelayed overtaking on lane i' = i+1 without any change of velocity 
{y' = v) by vehicles which would be hindered by slower vehicles on lane i. Analogously, term 
( |22b ) reflects undelayed overtaking maneuvers on lane i' = % — 1. Term ( [22cj ) delineates 



situations where a vehicle cannot be immediately overtaken by a faster vehicle so that the 
latter must stay on the same lane {i! = i) and decelerate to the velocity v' — w of the 
hindering vehicle. 

B. Lane-Changing 

We come now to the specification of the lane-changing term (df>i/dt)\ c . This has the form 
of a master equation: 

(fr) = E W 1 (t\ l ')p l ,(r,v,v ,t) (23a) 

\ / lc i'&i) 

- W 1 (i'\i)p i (r,v,v ,t). (23b) 
i'&i) 

Term (|23a|) describes an increase of phase-space density pi(r, v,v ,t) due to changes from lane 
i' i to lane i by vehicles with actual velocity v and desired velocity v . The frequency of 
lane-changing maneuvers is proportional to the phase-space density pii(r,v,vo,t) of vehicles 
which may be interested in lane- changing. Analogously, term (|23b| ) reflects changes from lane 
i to another lane %' causing a decrease of pi(r, v, v , t). For the corresponding rate Wi(z'|z) of 
lane- changing maneuvers we assume 

Wi(i'\i) = i[l -Ci(r,v,v ,t)}8i^i +1 ) (24a) 

i 

+ J- [1 - Q(r, f , f o, t)R'(i-i) (24b) 



Term fl24a| ) describes lane changes by vehicles which are able to escape a queue on the 



left-hand lane i' — i + 1. is the average waiting time until this is possible. Analogously, 
term ( |24b| ) reflects delayed overtaking maneuvers by impeded vehicles on the right-hand lane 
z' = z — 1. For free traffic flow on European autobahns we must set 

i = (25) 

i 

since overtaking is only allowed on the left-hand lane i' — % + 1. 



Term (24c) delineates spontaneous lane changes to lane £' = z + 1 by vehicles that desire 
to move on lane iq > i. Analogously, (|24d|) corresponds to changes to lane %' = i — 1 by 
vehicles that desire to move on lane io < i. and T[~ denote the related average waiting 
times. 

For American freeways we have 

Ti = n(r, v, v , t) = Ci(r, v, v , t) , (26) 

since lane changing by queuing vehicles corresponds to delayed overtaking which is described 
by terms ( |24a|) and ( |24b|) . For free traffic flow on European autobahns we have pf = due 
to a regulation prescribing to use the right-most lane if possible. This also applies to queuing 
vehicles so that, additionally, 

n = l . (27) 

We close with some remarks regarding the order of magnitude of the waiting times Tf~ 
and T^ 1 . First of all, they depend on the vehicle density on the neighboring lane % ± 1. 
Moreover, for the left-most lane we must set 

^T = and 4^ = 0, (28) 



whereas for the right-most lane we always have 



On American freeways, the waiting times Tf~ for spontaneous lane changing are greater than 
the waiting times for delayed overtaking since overtaking is usually carried out as soon 
as possible, whereas changing towards the desired lane is done at one's convenience (if one 
does not intend to leave the next exit), i.e. 

11 

l i 1 i±l 

On European autobahns, the waiting time T~ for spontaneous lane changing (to the right- 
hand lane) is even much larger than the waiting time T^i 1 for delayed overtaking (on the 
left-hand lane) since the right-hand lane is only accepted when it offers really large gaps. 
(Nobody wants to move on the slower lane.) 
Defining the abbreviations 

1 1_Ci + ^ (31) 



and 



± i J-i J-i 
we can rewrite expression ([24]) in the simple form 



n t+ t + 



1 1 — C; pfri 

+ , (32) 



W x {i'\i) = + ^r^i'(i-i) • (33) 



C. Distinction of different vehicle types 

The gas-kinetic traffic model developed above can be easily generalized to cases where 
effects of different vehicle types (e.g. cars, trucks) or driving styles (e.g. aggressive, timid) 
are to be investigated. Then, the phase-space density pi(r, v,vo,t) can be splitted into partial 
phase-space densities p1(r,v,v ,t) describing driver-vehicle units of type a: 

Pi(r,v,v ,t) = ^p 4 a (r,f,f ,t) . (34) 



The corresponding gas- kinetic equations have an analogous form to @ with (0), (|TTD, 
(0), and They read 

+ <9r lPi j + l Pt 4 T a ) 



EE J dv' J dw J dw W£ b (v, i\v', i'; w, i')p^(r, v ', t> , typ^r, w, w , t) 

b V w<v> 

^2YlJ dv 'J dw J dw W2 b (v',i'\v,i;w,i)p^(r,v,Vo,t)p b i (r 1 w,w ,t) 



b i w<v 



+ E W?(i\i')$(r,v,v ,t)-J2 W?(i'm(r,v,v ,t) 

+ 0? + (r, v, v , t) - z>r(r, v, v , t) . (35) 

Here, we have taken into account that the relaxation time T a and the reasonable distribution 
of desired velocities may depend on driving style or vehicle type a. The specification of 
the interaction rates and the lane-changing rates is analogous to the previous 
discussion (cf. (gD and (|D). 



III. DERIVATION OF MACROSCOPIC TRAFFIC EQUATIONS 

The gas-kinetic traffic equations are not very suitable for computer simulations since they 
contain too many variables. Moreover, the phase-space densities are very small quantities 
and, therefore, subject to considerable fluctuations so that a comparison with empirical data 
is difficult. However, the special value of gas-kinetic traffic equations is that they allow a 
systematic derivation of dynamic equations for the macroscopic (collective) quantities one 
is mainly interested in. 

A. Definition of Variables 

The most relevant macroscopic quantities are the vehicle densities 

Pi(r,t) = J dv J dv pi(r,v,v ,t) (36) 

and the average velocities 



Vi(r,t) = (v)i = J dvvPi(v;r,t) (37) 

on lanes i. Here, we have applied the notation 

F t (r,t) = (f(v,v ))i = fdvf dv f(v,v ) Pi{r, '" ,V ° ,t) (38) 

J J Pi{r,t) 

and introduced the distribution of actual velocities 

pf + \ f A M r > u > u o,f) Pi(r,v,t) . , 

on lane i. Analogous quantities can be defined for vehicles entering and leaving the road at 
entrances and exits respectively. 

"t( r , t) = Jdvj dv uf(r, v, v , t) (40) 
are the rates of entering and leaving vehicles, and 

V?(r,t) = (v)f = J dvvPt(v;r,t) (41) 
their average velocities, where 

are the velocity distributions of entering and leaving vehicles respectively. In addition, we 
will need the velocity variance 

e^t) = ([v - = Jdv[v- V t (r,t)} 2 P t (v;r,t) = (v 2 ), - ((v),) 2 (43) 

and the average desired velocity 

V 0i (r,t) = [dv [ dv v M r '"' V Z' t) (44) 
J J pi{r,t) 

on each lane % as well as the average interaction rate 

11/". r 

dv pi(r,v,t) dw (v — w)pi(r,w,t) (45) 



T? Pi (r,t) 

of a vehicle on lane % with other vehicles on the same lane. 

Finally, we will assume that the proportion q of freely moving vehicles is a function of 
the density pi on the respective lane i: 



Ci(r, v, v , t) = Ci[pi(r, t)] . (46) 

The probabilities with which slower vehicles can immediately be overtaken will mainly de- 
pend on the density on the neighboring lane i ± 1, i.e. 

Pt = Pt{pi±i) • (47) 
A similar thing holds for the waiting times: 

ft = ft( Pl ±i) , ft = ft{Pi±i) ■ (48) 

Consequently, 

pf = pf(j>i-i, Pi, Pi+r, r > /) (49) 

and 



Tt = Tt( Pi , Pi±1 ;rJ). (50) 

From a theoretical standpoint, a generalization to more complicated functional dependences 
is easily possible, but their determination from empirical data becomes very difficult, then. 

B. Derivation of Moment Equations 

We are now ready for deriving the desired macroscopic traffic equations from the gas- 
kinetic equation (|) with (0), ®, (HI]), ©, (||), and (||). Integration with respect 



to Vq gives us the reduced gas-kinetic traffic equation 



- (1 - Pi)p~i(r, v, t) J dw(v - w)p~i(r, w, t) (51b) 
+ pt 1 pi-i{r,v,t) dw(v-w)p i „ 1 (r,w,t) (51c) 



+ Pi+iPi+i( r , v ,t) J dw (v - w)p i+ i(r,w,t) (51d) 
- (pf + Pi )Pi(r, v, t) J dw(v - w)p~i(r, w, t) (51e) 



1 ~ / \ 1 _ , 

+ 7f ^p t -i{r,v,t) - ■=+pi{r,v,t) 
1 ~ , 1 -/ 

+ 7f—pt+i{r,v,t) - —pi(r,v,t) 
+ vt(r,v,t) - uf(r,v,t) 



(51f) 

(51g) 
(51h) 



with 



and 



Voi{v) = V 0i (v;r,t) = / dv v 



Pi(r,v,v ,t) 
Pi(r,v,t) 



t) = J dv vf(r,v,v ,t) . 



(52) 



(53) 



In formula (|5l"D, the deceleration term ( |51b|) stems from ( p2cj ), the terms (|51c|) to ( |51e| ) 
reflecting immeditate overtaking come from (|22a|) and ( |22b|) , and the lane-changing terms 
(15111) , (|51g|) originate from (|33]). The adaptation term (dpi/dt)^ yields no contribution. 

We will now derive equations for the moments (v k ) by multiplying (|51|) with v k and 
integrating with respect to v. Due to 



i k d (~ 
dv\ PiC ' 



Voijv) ~ v S 
T 



dv kv k 1 I pic 



Voijv) ~ ^ 
■; T 



kci 



(54) 



dvv 



dv 2 



i— J = k(k - l)*-jr{v )i 



(55) 



and 



(1-Pi) dvpi(r,v,t) dw(wv k -v k+1 )p l (r,w,t) = (l-p i )( Pl ) 2 ((v) l (v k ) t ~(v k+1 ) l ) (56) 



we obtain the macroscopic moment equations 



(v k ~ 1 vo)i - (v h )i + (k — l) — (v k ) 



Ci 



+ (i-p i )(p i n(v) i (v k ) i -{v k+i ) i ) 



+ ^Pi-i[(v k )i-i + Kti) - %Pi[{v k )i + K*) 



Pi 



1 i+1 1 i 

+ ^(r,t)(v k )t-^(r,t)(v k )-, (57) 



[v k )f = Idvl dv v k ^^ = (dvv k ^4 (58) 



where 



y\ (r,t) J v i (r,t) 

and 

T° r r 
fci = — dvpi(r,v,t) dwv k (v - w)p~i(r,w,t) - (v k )i. (59) 
Pi J J 



C. Fluid-Dynamic Multi-Lane Traffic Equations 

In order to derive dynamic equations for the densities pi and average velocities Vi, we 
need the relations 

(v\ = ([Vi + (v- Vi)] 2 ), = (V{) 2 + 2V t (v - Vi)i + ((v - V) 2 ) t = {Vf + Bi (60) 

and 

Pi ~dt ■ = m^ {v) - ] - Vi m ■ (61) 

Applying these and and using the abbreviations 

4 = f| + ^E, (62) 

' i i i 

equation (|57|) gives us the density equations 

^F + V ^ = - p ^t + ^ (r ' t} " ^ (r ' £) (63a) 

+ Pl_l_Pl_ + P2±l_PL. (63b) 

T i-1 T i T i+l T i 

This result is similar to previous multi-lane models. In addition, we obtain the velocity 
equations 



dV. dVi 



dt 



Or 



d{pA) + C -^(V 0i -V i )-(l- Pi )( Pi )% 



dr 

Ti-l 



T 

v: 



pi+i 



(V i+1 - Vi) 



i Pi- 1 n V 1 i P1+ 1 „ V 1 
+ T^~Pi- lh -i-l + Tf^—Pi+l^i+l 

+ 4(V+-V i )-vi(Vr-V i ) 



Pi + Pi 



(64a) 
(64b) 

(64c) 
(64d) 



after some lengthy but straightforward calculations. The terms containing the rates vf 
and v~ reflect entering and leaving vehicles, respectively. Whereas the terms ( |63a| ) and 
( |64a|) correspond to the effects of vehicle motion, of acceleration towards the drivers' desired 
velocities, and of deceleration due to interactions, the terms ( 63b ), (|64b|) , and (|64cj ) arise from 
overtaking and lane-changing maneuvers. ( |64b ) comes from differences between the average 



velocities on neighboring lanes. In contrast, ( |64cj ) originates from the fact that overtaking 
vehicles are, on average, somewhat faster than the vehicles passed. Whereas (|64c|) tends to 
produce differences between the average velocities of neighboring lanes, these are reduced by 
( |64b| ) . The term ( 64d|) has a similar form and interpretation like (|64b ). It is only negligible if 
entering vehicles are able to adapt to the velocities on the merging lane and exiting vehicles 
initially have an average velocity similar to that on the lane which they are leaving so that 

In order to close equations (|63]) and fl64j) , we must specify the interaction rates l/T^ , the 
correction terms JCj, and the variances 9{. Utilizing that the empirical velocity distributions 
Pi(v;r,t) are approximately normally distributed \ fT] , 25 ,^9 j |52| , we have 



Pi(v;r,t) 



1 c -k-K;(r,f)] 2 /[28,;(r,t)l 



(65) 



which implies 



(66) 



and 



0; 



(67) 



With a detailled theoretical and empirical analysis it can be shown [2C] that the variances 
8i can be well approximated by equilibrium relations of the form 



u i = r~{ — j7 — \ > ( bi= v 

where Ci(pi,Vi) = ((v—Vi)(vo—Voi))^ denotes the density- and velocity-dependent equilibrium 
covariance between actual velocities v and desired velocities vq on lane i. In the case of a 
speed limit there is v V so that « 0. 

For the average desired velocities we have 

^ = y O i(r,0~^(r,/) (69) 

since 

Ai(«o;r,/) -PoiKM) ^0 (70) 

due to the smallness of T T . 

IV. EFFECTIVE MACROSCOPIC EQUATIONS 

We will now investigate how the effective macroscopic traffic equations for the total 

cross-section of the road look like. For this purpose we introduce the effective (average) 
density 

p(r,t) = jJ2Pi(r,t), (71) 

i 

where / = I(r,/t) is the number of lanes at place r, and calculate the average Fj over all 
lanes for every lane-specific quantity Fj according to 

F = Fi-Hj- p Fi - (72) 



This implies 



pF = jJ2^- (73) 



In addition we define the effective rates 

^M) = yE^M) (74) 

i 

of entering and leaving vehicles, respectively, and the corresponding lane-averages 



F± = Ft = Y,^Ft- (75) 



In the following, we will apply the factorization approximation 



Pi F i G i K Pi F i G i = pFG. (76) 



Summation of equation (|57|) over i and division by / gives us 
d , / i \ \ d , , r— tt N kc 



oibWh) + fr(p( vk+1 h) = ^[(^Mi - <«*>« + ( fc - i)(A/c,) («*>«] 



+ ^ + (r,t)(^)+-^(r,t)(^)-, (77) 
since the lane-changing terms cancel out each other: 

7 E 4-«n<«*>*Fi " 7 E ;4*<A = 7 E 4p 4 -<A- - 7 E :W fc >* = • (78) 

A similar thing holds for the terms containing the corrections K.f. As a consequence, only the 
terms describing vehicle motion, acceleration towards the drivers' desired velocities, velocity 
fluctuations, deceleration due to vehicle interactions, and entering or leaving vehicles are 
remaining. 

Now, apart from the lane-average 9 = 9i of the velocity variances 9i we additionally 
define the total effective variance 



& = ((v- vy), = ([(v - vs) + (Vi - v)Y) % = e + {Yi- vy . (79) 

Applying this and relation 



(v 2 ) i = (V i ) 2 + 9 i =[V+(V i -V)] 2 + 9 



= V+ (Vi - V) 2 + 9 = V +Q (80) 

to equation ([77]), after some lengthy but straightforward calculations we obtain the effective 
density equation 

dt dr P dr 
and the effective velocity equation 

8V BV d(pQ) C P/jr t r\ m \ 2n 

+ u +(y + -V)-v-(V~ -V). (82) 



A. Comparison with former macroscopic traffic models 



Taking a look at the effective density equation (BTJ), we recognize that it has the form 
of a continuity equation. This reflects the conservation of the number of vehicles, i.e. no 
vehicle on the considered road is produced or gets lost. The continuity equation was first 
proposed in the fluid-dynamic traffic model by Lighthill and Whitham P|pT|, and it is also 
a component of most other macroscopic traffic models |],^-[TT|. 

The effective velocity equation (|8"2|) can alternatively be represented in the form 



pd ~L + pV W = ~W + r {Ve ~ V) + V+{V+ ~ V) - V ~ {V ~ - V) ' (83) 



where we have introduced the so-called traffic pressure [^q , ^ ,p3 



v = P e. 



(84) 



The effective relaxation time 



r{ P ) 



T 



(85) 



increases with decreasing effective proportion c(p) of freely moving vehicles and reflects 
the well-known "frustration effect" that the average relaxation time increases with growing 
vehicle density p. In addition, we have introduced the equilibrium velocity 



V e = V -r(l-p)p6 



(86) 



which is determined by the average desired velocity Vq, diminished by a term arising from 
deceleration maneuvers due to vehicle interactions. The velocity equations proposed by Pri- 
gogine et al. [46|, Paveri- Font ana |47j], Phillips |]52| , |53| , Payne ||, as well as Kerner and 
Konhauser can all be written in the form (|83| ) |TT| . However, the specification of the 
functions V, r, and V e is varying from one model to another. The terms i/ ± (V A± — V) have 
been neglected by all previous models, despite the fact that entrances (and exits) reduce the 
freeway capacity due to V ± < V . Moreover, the effective gas-kinetic models proposed by 
Prigogine et al. jPJ, Paveri- Font ana fi"? |, and Phillips |52| , |53| imply V c = Vq — r(l — p)pQ 
instead of relation (|36"D. This is a consequence of the fact that they treat the total freeway 



section like a single lane of greater capacity and possibilities for overtaking. However, the 
approximation « 9 is only valid for V, ~ V. Although it is normally justified for American 
freeways, on European autobahns it is only valid at densities greater than 35 vehicles per 
kilometer and lane |20| . 



V. SUMMARY AND OUTLOOK 

In this paper we have derived a macroscopic traffic model for uni-directional multi-lane 
roads. Our considerations started from plausible assumptions about the behavior of driver- 
vehicle units regarding acceleration, deceleration, velocity fluctuations, overtaking, and lane- 
changing maneuvers. In addition, we distinguished freely moving and queuing vehicles since 
these behave differently. The resulting gas-kinetic traffic model is a generalization of Paveri- 
Fontana's Boltzmann-like traffic equation. It can be extended to situations where different 
vehicle types or driving styles are to be investigated. 

The gas-kinetic traffic equations not only allow to derive dynamic equations for the 
vehicle density on each lane, but also for the average velocity. In addition, we obtained 
effective macroscopic traffic equations for the total cross-section of the road from the multi- 
lane model. The resulting equations implied corrections with respect to previous traffic 
models, but agree with the ones by Prigogine et al. and Paveri-Fontana in the special case 
Vi » V « V ± . 

Strictly speaking, the presented model is only valid for small vehicle densities. However, 
with the methods discussed in Refs. fllj20l it can be easily generalized to a traffic model for 



arbitrary densities. Since the calculations are rather lengthy but straightforward, they have 
been omitted, here. 

Present research focuses on the implementation of the multi-lane model on a computer 
with the objective of carrying out numerical traffic simulations. The functional relations re- 
garding the probabilities of overtaking, the proportions of freely moving vehicles, the waiting 
times for lane-changing maneuvers, etc. are evaluated on the basis of empirical data. Some 
of the questions we are going to investigate by means of multi-lane simulations are the 
following: 



1. In which way does on-ramp traffic influence and destabilize the traffic flow on the other 
lanes? How does the destabilization effect depend on the traffic volume, the length of 
the on-ramp lane, the total lane number, etc.? 

2. In case of a reduction of the number of lanes, is it better to close the left-most or the 
right-most lane? 

3. Is the organization of American freeways or of European autobahns more efficient, or 
is it a suitable mixture of both? Remember that American freeways are characterized 
by uniform speed limits and the fact that overtaking as well as lane changing is allowed 
on both neighboring lanes. In contrast, on European autobahns often no speed limit is 
prescribed (at least in Germany) and average velocity normally increases with growing 
lane number since overtaking is only allowed on the left-hand lane. 

4. In which traffic situations do stay-in-lane recommendations increase the efficiency of 
roads? 
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